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We report on a study of (001) oriented fee Fe-Ni alloy surfaces which combines first-principles 
calculations and low-temperature STM experiments. Density functional theory calculations show 
that Fe-Ni alloy surfaces are buckled with the Fe atoms slightly shifted outwards and the Ni atoms 
inwards. This is consistent with the observation that the atoms in the surface layer can be chemi- 
cally distinguished in the STM image: brighter spots (corrugation maxima with increased apparent 
height) indicate iron atoms, darker ones nickel atoms. This chemical contrast reveals a c(2 x 2) chem- 
ical order (50% Fe) with frequent Fe-rich defects on the Fe64Ni3e(001) surface. The calculations also 
indicate that subsurface composition fluctuations may additionally modulate the apparent height of 
the surface atoms. The STM images show that this effect is pronounced compared to the surfaces of 
other disordered alloys, which suggests that some chemical order and corresponding concentration 
fluctuations exist also in the subsurface layers of Invar alloy. In addition, detailed electronic struc- 
ture calculations allow us to identify the nature of a distinct peak below the Fermi level observed 
in the tunneling spectra. This peak corresponds to a surface resonance band which is particularly 
pronounced in iron-rich surface regions and provides a second type of chemical contrast with less 
spatial resolution but one that is essentially independent of the subsurface composition. 

PACS numbers: 73.20.-r, 68.37.Ef, 73.61.At 

Keywords: nickel; iron; alloy; surface electronic structure; scanning tunneling microscopy; density functional 
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I. INTRODUCTION 

The properties of alloy surfaces are of great impor- 
tance for surface- and interface-related phenomena like 
segregation and catalysis, furthermore, examination of 
the surface makes it possible to study bulk properties not 
easily accessible by other methods. Scanning tunneling 
microscopy (STM) and spectroscopy (STS) have become 
indispensable methods for imaging alloy surface*^^ 
and surface alloys^i 7 - with atomic resolution and chemi- 
cal contrast. This paper focuses on the interpretation of 
STM/STS data of a transition metal alloy, the Fe64Ni36 
alloy, based on first principles calculations. In the course 
of this analysis two separate modes of chemical contrast 
in STM images are discussed: First, apparent height dif- 
ferences of different chemical species providing chemical 
contrast in an atom-by-atom fashion; Second, chemical 
contrast based on specific features of the tunneling spec- 
tra, in particular surface state peaks. 

On many transition metal surfaces, surface states and 
resonances near the Fermi level have been detected by 
STS-—————. However, according to experience and 
the most basic theory of the STMi^ only a small fraction 



of surface states, specifically those which are sufficiently 
delocalized perpendicular to the surface, can by seen by 
the STM, i.e., contribute to the tunneling current. A 
class of surface resonances accessible by STS is that with 
s — p z — d z 2 symmetry at the T-point observed on the 
bcc(001) surfaces of Fe^, Cr— , V—, W— , and their sur- 
face alloys Cr-Feii and Si-Fe^, and on the hcp(0001) 
surface of Gd/W(110)i£. These states may serve as an 
extra source of chemical contrast as alloy components 
tend to locally quench the surface resonance as demon- 
strated for the dilute surface alloy Cr/Fe(001)ii and the 
chemically ordered surface alloy Si/Feg6.4Si3. e(Q01)^. 

In this article we apply these ideas to the fee alloy sur- 
face Fe64Ni36(001) which is different from the mentioned 
surface alloys in at least two ways: (1) Subsurface com- 
position fluctuation a 17 ' 18 lead to relatively strong varia- 
tions of the apparent height of the surface atoms on the 
nanometer scale, which makes the chemical assignment 
of the individual atoms, successfully done for the alloy 
surfaces Pt 2 5Ni 75 (lll)i or Pt 5 oRh 5 rj(001)i, rather diffi- 
cult. (2) This surface does not show clusters or chains of 
Fe or Ni atoms which a priori could be considered suf- 
ficiently large to support ID or 2D surface resonances 



characteristic for the elements constituting the alloy. 

Fe64Ni 36 has a very low thermal expansion for tem- 
peratures up to its Curie temperature of approximately 
500 K and is therefore known under the trademark Invar 
since its discovery more than 100 years ago. Invar and 
other iron-rich Fe-Ni alloys are applied in numerous ther- 
mally stabilized or matched components. The thermal 
anomalies of these materials arc explained by a tempera- 
ture dependent atomic volume of Fe in the ferromagnetic 
state counteracting "regular" thermal expansion due to 
lattice anharmonicity (cf., e.g., Refs. [T7H19II20I ) . A non- 
collinear magnetic ordering has been predicted^ 7 , and ob- 
served recently 2 ^, which depends on composition fluctu- 
ations and chemical orde r 17 ' 18 ! 22 . Therefore, local prop- 
erties like chemical order, surface and subsurface compo- 
sition, and surface electronic structure are very relevant 
to deepen the understanding of this technologically im- 
portant system. 

After presenting the experimental results, we first in- 
vestigate the surface geometry of Fe-Ni(OOl) alloys from 
first-principles, in particular, we will demonstrate that 
the alloy surface is buckled and this buckling allows us 
to explain the chemical contrast observed in atomically 
resolved STM images. Second, we present a detailed first- 
principles investigation of the surface electronic structure 
of the Fe64Ni36(001) alloy surface and show that the sur- 
face resonance found just below the Fermi level of Fe-rich 
surfaces is in accordance with the local spectroscopic fea- 
tures found in the STM experiment. 
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II. EXPERIMENT 



A. Setup and sample preparation 

The experiments were done using a low-temperature 
scanning tunneling microscope (LT-STM) mounted in- 
side of cryo-shields cooled by a bath cryostat filled with 
liquid N 2 or liquid He (Omicron LT-STM). The micro- 
scope is operated with electrochemically etched W tips 
conditioned by sputtering with Ar + under ultrahigh- 
vacuum (UHV) conditions. The LT-STM system is 
equipped with an Auger electron spectroscopy (AES) sys- 
tem for a surface sensitive chemical analysis. A separate 
UHV chamber was used for low energy electron diffrac- 
tion (LEED) measurements at 300 K and above to probe 
the temperature dependence of the surface chemical or- 
der. The pressure during measurement in the separate 
STM chamber of the LT-STM system is around 1 x lO" 11 
mbar, and in the chamber with the LEED system in the 
10~ n mbar range. 

The Fe64Ni36(001) single crystal was typically pre- 
pared by sputtering with 2 kcV Ar + ions and subsequent 
annealing at 770 K for 10 minutes. 



FIG. 1: (a) STM image of the Fe 6 4Ni 36 (001) surface (sam- 
ple voltage -2 mV, tunneling current 8 nA) . The three details 
highlight a c(2 x 2)-ordered area (unit cell marked by square) 
and both a short Ni-rich and a short Fe-rich anti-phase do- 
main boundary segment, (b) Profile of the STM image along 
the white line in (a). The dashed average curve indicates 
a long-wavelength buckling on the nanometer scale in con- 
trast to the short-wavelength buckling distinguishing Fe and 
Ni atoms, (c) Map of automatically detected Fe atom loca- 
tions (corrugation maxima) in the STM image omitting the 
Ni atoms, with each of the Fe atoms assigned to one of the 
two (2 x 2) anti-phase sublattices (open and full circles, occa- 
sional oval shapes are due to the detection of more than one 
corrugation maximum per atom caused by image errors). 



B. Chemical contrast and surface order 

Figure 1 shows an STM image of the Fe 64 Ni 36 (001) 
surface at 80 K. The distinctly different brightness (ap- 
parent height) of adjacent atoms indicates a chemical 
contrast between Fe and Ni atoms. Images from similar 
Fe64Ni36 surfaces covered by additionally deposited Fe 
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(about 1 monolayer) in order to artificially increase the 
Fe content suggest that the "dark" species is Ni. This 
atom-by-atom chemical contrast reveals that the surface 
is locally ordered with a c(2 x 2) superstructure [top right 
detail in Fig. 1(a)]. Since the Fe surface concentration is 
higher than 50% (cf. next paragraph), the excess Fe has 
to be inserted as Fe-rich defects, which often are short 
pieces of Fe-rich anti-phase domain boundaries (bottom 
left detail). Occasionally, however, also Ni-rich defects 
are visible in the STM images (middle detail) , indicating 
that the disorder is larger than required to accommodate 
the excess Fe. A LEED experiment probing the tem- 
perature dependence of the (1/2, l/2)-type superstruc- 
ture spots normalized to (1, 0) integer spots shows a re- 
versible rather abrupt disappearance and reappearance of 
the signal near 550 K indicating an order-disorder tran- 
sition there. Atomically resolved STM images with up to 
20 x 20 nm frame size suggest typical anti-phase domain 
diameters between 1 and 10 nm at 80 K. 

The chemical analysis by automatically counting 
bright and dark atoms based on their apparent height is 
very unprecise (46%±12% Ni) since the histogram of ap- 
parent atom height values does not show a clear bimodal 
distribution of dark and bright atoms. The reason is a 
significant buckling of the surface in the entire spatial 
wavelength range of the image [dashed line in Fig. 1(b)], 
interfering with the buckling providing the chemical con- 
trast between adjacent atoms. Only strong high pass 
filtering (not shown in Fig. 1), leaving only height differ- 
ences of adjacent atoms in the image, leads to separable 
peaks in the histogram of apparent atom heights, which 
can be fitted by a double gaussian. Combining the re- 
sulting Ni concentration with results from manual atom 
counting and a quantitative AES analysis, using the sur- 
face sensitive low-energy peaks of Fe and Ni, results in a 
Ni concentration of 41% with an estimated error of ±2% 
for the annealed Fe64Ni36 surface — a small Ni enrichment 
relative to the nominal bulk concentration of 36%. 

The buckling of the surface layer on the nanometer 
scale may be explained by the subsurface distribution of 
larger Fe and smaller Ni atoms. Since there exists no 
long-range order in the Fe64Ni36 bulk, the first hypoth- 
esis assumes a fully random order of the subsurface Fe 
and Ni atoms, which implies random composition fluc- 
tuations. However, the surfaces of the disordered alloy 
Pt25Ni75(lll)i, which shows weak surface order, or of the 
alloy Pt25Rh 7 5 (100)2 do not show a comparable surface 
buckling on the nanometer scale. Therefore most likely 
the second hypothesis applies which assumes subsurface 
local order as observed in the Fe64Ni36 bulk by element 
specific x-ray scattering experiments^. This would be 
an arrangement of subsurface Fe and Ni atoms which is 
rather similar to that seen on the surface with c(2 x 2) 
ordered regions (50% Fe) interrupted by frequent small 
Fe-rich domain boundary fragments. The corresponding 
composition fluctuations would then have a characteristic 
length scale above 1 nm, i.e., the diameter of Fe-rich de- 
fects or ordered domains, as observed, while fluctuations 



in randomly distributed alloys are strong only at very lo- 
cal scales, which have less influence on the buckling of ad- 
jacent atomic layers. We remark, however, that it is likely 
that the chemical order on the surface is higher than 
that in the subsurface layers due to smaller barriers for 
atom exchange on the surface and the low order-disorder 
transition temperature of the system around 550 K, i.e., 
a temperature which seems prohibitive for efficient bulk 
diffusion (assuming equal transition temperatures for the 
surface and the bulk). 

As mentioned at the beginning of this section, this 
convolution of surface and subsurface effects make the 
chemical analysis of the surface based on apparent atom 
heights rather cumbersome. In addition, in many cases 
the STM images do not show the individual metal atoms 
and the brightness variations cannot be unambiguously 
assigned to either surface chemical contrast or surface 
buckling. However, tunneling spectroscopy may be used 
to separate surface and subsurface effects - although with 
less spatial resolution and less precision as demonstrated 
in the next section. 



C. Chemical contrast by scanning tunneling 
spectroscopy 

Figure 2 shows STS images and tunneling spectra of 
the Fe64Ni 36 (001) surface at 5 K and spectra only for a 
surface at 80 K. Data acquisition was done by measuring 
I/V curves for every image point with the tip position 
zq adjusted to 1 nA at -1 V tunneling voltage in every 
image point [Fig. 2(a)], followed by numerical differenti- 
ation. The spectra shown in Fig. 2(c) are averages over 
a few tens of individual spectra from equivalent image 
points of the same STS image without additional smooth- 
ing. Both experiments, 5 and 80 K, show a distinct peak 
at 0.35±0.1 eV below the Fermi edge, which is intense, 
however, only in a very small fraction of the surface. The 
first principles calculations in the following sections will 
describe the character of the underlying electronic states 
in detail. In this section we will use the term "surface 
resonance peak" for the sake of simplicity. The small 
peak shift of 0.1 eV between 5 and 80 K is not significant 
within the precision of this measurement. The small peak 
visible in the 80 K measurement [arrow in Fig. 2(a)], ap- 
pears in the spectrum of every image point in this data 
set only, which is characteristic for a tip-related state. 

The peak shape is rather symmetric, indicating lit- 
tle dispersion of that part of the surface resonance band 
that is visible for the STM (cf. Sect.[mB|. The width 
(FWHM) is about 250 meV, which is quite large com- 
pared to similar peaks on the Cr(001)^ or Gd(0001)^ 
surface at 5 K. However, it is comparable to the width 
of the resonance of single Fe atoms on the Pt(lll) sur- 
face at 5 K (rj 300 meV)^. The calculations in the next 
sections show that the peak can indeed be assigned to 
Fe-rich defects on the chemically ordered Fe64Ni36(001) 
surface. 
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height of the surface does not obviously correlate with 
the map of the surface resonance intensity. Although a 
weak correlation should exist, because the composition of 
the surface monolayer contributes to the apparent height 
as well, it is too weak to be relevant compared to the 
height variation caused by the compounded effect of con- 
centration inhomogeneities in each of several monolayers 
below the surface monolayer. Thus, the tunneling spec- 
tra provide a method to map the chemical information 
independent of the geometric height information on ev- 
ery image scale without the requirement to resolve the 
surface atomically. We remark that, although desirable, 
it would be extremely difficult to combine such a spec- 
troscopic analysis with the simultaneous acquisition of an 
atomically resolved image, mainly due to the instability 
and short lifetime of very sharp tips. 



III. DFT CALCULATIONS 



A. Surface relaxation 



FIG. 2: (a) Constant current image of the Fe64Ni36(001) sur- 
face at 5 K, 1 nA tunneling current, and - IV sample voltage, 
(b) Corresponding current map at -0.4 V, showing the in- 
tensity of a surface resonance below the Fermi-level (shaded 
area in 5 K spectrum to the right), (c) dI/dV(V) spectra at 
5 K corresponding to the images at left and of a separately 
prepared sample at 80 K. Shown are spectra near extremal 
points of the surface resonance intensity, both for maximum 
(A) and minimum (B) intensity. The small peak at -0.1 V in 
the 80 K spectra is a tip-related state (arrow). The contour 
lines in (a) and (b) encompass the relatively small surface 
fraction that shows spectra which are more similar (smaller 
root mean square deviation) to spectrum A than to spectrum 
B. The remaining larger surface fraction shows spectra which 
are more similar to spectrum B. Evidently, the surface frac- 
tion with significant surface resonance intensity is rather small 
(26% in this data set) consistent with the assumption that it 
marks Fe-rich regions. 



Mapping the tunneling current at a voltage just below 
the low-energy edge of the resonance peak by plotting 
I(x, y, Zo,— 0.4 V) from the complete I/V dataset pro- 
vides a map of the surface resonance intensity acquired 
in a distance of more than 0.5 nm from the surface plane 
[Fig. 2(b)]. By evaluating the similarity of the spectrum 
in each image point to the peaked spectrum (A) and the 
"flat" spectrum (B) wc find (see caption Fig. 2) that the 
surface resonance peak is a minority property of the sur- 
face, with the larger surface fraction showing no or only 
a small peak, suggesting a correlation to the Fe-rich de- 
fects seen in the atomically resolved images [Fig. 1(a)]. 
Direct comparison of the classification by the tunneling 
spectra [contours in Fig. 2(a,b)] with the constant cur- 
rent image [Fig. 2(a)] demonstrates that the apparent 



We first investigate the surface geometry of the Fe-Ni 
alloy. The relation between crystal structure and mag- 
netic ordering is of crucial importance for iron-based al- 
loys. The calculations have been done using the full- 
potential linearized augmented plane wave method^ 
(FP-LAPW) code FLAIRS. Two approximations for 
the density functional potential have been used: the local 
density approximation (LDA) according to Vosko, Wilk, 
and Nusaii— and the generalized gradient approximation 
(GGA) according to Perdew, Burke and M. Ernzerho^. 
Spherical wave and plane wave expansions were trun- 
cated at Imax = 8 and E cut = 11.6 Ry respectively. The 
difference between the last two charge density iterations 
was less than 10 -6 e a.u.~ 3 . In all the surface geometry 
optimizations, atomic positions in the surface and two 
subsurface layers were fully relaxed^. The final relaxed 
structure was assumed when the force on each atom was 
smaller than 5 meV/A. Nine-layer ordered films have 
been used to simulate fcc(001) alloy surfaces with differ- 
ent ordering. For all the calculations we have used the 
experimental lattice constant a — 3.58 A2£. Minimizing 
the total energy as a function of lattice constant for a 
NiFe3 crystal yields the value a = 3.60 A. Our model 
geometries are derived from a NiFe3(001) surface struc- 
ture composed of c(2 x 2) NiFe(001)-layers alternating 
with pure Fe(001) layers. It should be noted here that 
this corresponds to the CU3AU structure, also known as 
LI2. Our models are in accordance with STM experi- 
ments, which indicate approximately the same amount 
of iron and nickel atoms in the top surface layer. All 
model systems exhibit inversion symmetry and include 
18 atoms in the unit cell. Some models have only iron 
atoms in the subsurface layer and differ in the geometry 
of deeper layers, e.g. by interchanged Ni and Fe atom 
sites. We considered also systems where the subsurface 
layer includes both type of atoms. 
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Results of first principles atomic force minimizations 
for different models^ can be summarized as follows: 
the top surface NiFe ordered monolayer is buckled with 
Fe atoms pushed outwards and Ni atoms pushed in- 
wards. The buckling height in the top surface layer, 
0.07 to O.f 2 A, depends on the type of atoms in subsur- 
face layers. We may assume that also in the disordered 
fcc(OOl) alloy surface the Fe atoms are pushed outwards 
and their positions are determined by the local chemi- 
cal composition of subsurface atoms. The atomically re- 
solved STM images essentially reflect the core positions 
of the Fe-Ni (001) surface. The brighter spots can be 
associated with Fe atoms, the darker with Ni atoms. We 
conclude that the subsurface composition influences the 
top surface layer buckling as well as the first interlayer 
distance d 12 [Ad 12 = d u - d b f lk G (-0.07, +0.02) A}22. 
The presence of Ni in the subsurface layer results in a de- 
crease of the first interlayer distance and also the buck- 
ling of the surface is smaller. Different local chemical 
composition in the subsurface explains observed brighter 
and darker nanometer-sized areas in the measured STM 
images. 

Due to their reduced coordination, surface atoms in 
magnetic systems have higher magnetic moments as com- 
pared to the bulk^i. ft is well known that the surface 
magnetism depends not only on the coordination num- 
ber, but also on the details of atomic arrangement. For 
our models, we have found that the magnetic moments on 
the surface, 2.97 hb for Fe and 0.74 for Ni, are larger 
than the bulk values 2.61 \ib and 0.68 \xb- The given 
values arc local magnetic moments obtained by integrat- 
ing the charge distributions over the atomic (muffin-tin) 
sphere. The calculated work function for the alloy is 
4.50 eV. 



B. Surface electronic structure 

The peak in the tunneling conductivity found about 
0.35 eV below the Fermi level seems to indicate that a 
surface state or a surface resonance exists on the Fe64Ni36 
fcc(001) surface at this energy. In order to investigate the 
surface electronic structure, the spin-polarized local den- 
sity of states (LDOS) has been calculated for different 
surface geometries with FP-LAPW. We found a surface 
resonance below the Fermi level, but its energy and in- 
tensity was strongly dependent on the detailed film com- 
position. We needed thick film models (more than 30 
layers) to achieve well-converged results, rendering the 
computational demand for FP-LAPW too high. 

Therefore we have applied the tight-binding linear 
muffin-tin orbital method (TB-LMTO) in the atomic 
sphere approximation (ASA) to evaluate the surface al- 
loy features. The TB-LMTO code is based on the sur- 
face Green function formalism and incorporates realistic 
boundary conditions. The exchange-correlation poten- 
tial due to Vosko, Wilk, and Nusair— has been employed. 
The dipolc barrier effects due to charges spilling out from 



the sample into the vacuum are included in the formal- 
ism. The effect of randomness in the sample is described 
in the framework of the CPA (coherent potential approx- 
imation), which is able to treat arbitrary concentration 
profiles close to the sample surface. For further details 
see Refs. HH and HE 

The sample consists of two semi-infinite regions, the 
vacuum and the substrate, sandwiching the intermedi- 
ate region which also includes few 'substrate' and 'vac- 
uum' layers. The potentials in the intermediate region 
are determined selfconsistently while for the semi-infinite 
substrate bulk potentials are used. The vacuum region 
is represented by a flat potential; and the substrate is 
a binary alloy Fe64Ni36. The intermediate region con- 
sists of 10 atomic layers, namely of three empty-sphere 
layers adjacent to the vacuum region, two surface al- 
loy layers in which the composition can differ from that 
in bulk alloy, and five alloy layers in which we assume 
Invar bulk composition. In this way, charge densities 
in the intermediate region smoothly match those in the 
vacuum and substrate regions. The so-called principal 
layer o 32 ' 33 consisted of two atomic layers. The notation 
Fe a; Nii_ 2 ;/Fe y Nii_.y/Feo.64Nio.36 will be used to describe 
the composition of the two surface layers. We assume an 
ideal fee lattice with a lattice constant of a = 3.58 A for 
all layers. Thus, no lattice or layer relaxations are as- 
sumed. Changes of the interatomic distances due to the 
surface relaxation are only within few percent and do 
not affect the electronic structure of the Fe64Ni36 alloy 
substantially. The buckling of Fe and Ni atoms is also 
neglected. Possible errors introduced by this approxima- 
tion become small for Fe-rich alloy surfaces which are of 
interest here. 

The largest contribution to the tunneling current be- 
tween the tip and the sample comes from the regions 
around the point fcj| = T of the surface Brillouin zone, 
where the corresponding fcy -resolved local density of elec- 
tron states at given energy tends to decay most slowly 
into the vacuum (cf. eg. Ref. Q). Consequently, we 
start to analyze the spin- and k\\ -resolved local density 
of states, also called the Bloch spectral function (BSF), 
evaluated at the tip apex position for ku = T. 

The BSF for a pure iron surface layer on the Fe64Ni36 
alloy is shown in Fig. [3] The results clearly demonstrate 
that a surface resonance appears on iron-rich surfaces 
in an energy range centered at -0.3 eV, in good agree- 
ment with the experimental data. This surface reso- 
nance is formed by minority-spin electron states in the 
central part of the surface Brillouin zone and it has pre- 
dominantly d z 2 character. The totally symmetric Ai- 
component of the BSF is also shown. This Ai-component 
can be as well expressed as a sum of the s, p z , and d z 2 
orbital-resolved components. It is this particular compo- 
nent that gives rise to the surface resonance. The totally 
symmetric component is expected to decay into the vac- 
uum more slowly than the other components at the same 
energy do, and this behavior is indeed observed. Let us 
note that the Ai-component remains nonzero in the same 
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FIG. 3: Solid lines: Layer- and spin-resolved Bloch spectral 
functions at f for Fe/Feo.64Nio.36 (a) in vacuum 3.6 A above 
the surface, (b) in the top surface layer, and (c) in a layer deep 
in the bulk. Crosses denote the v4i-symmetry component of 
the BSF. 



energy range also inside the bulk substrate, as shown in 
Fig. [3jc) . The resonant electron states on the surface can 
thus hybridize with corresponding bulk states. 



(a) 



(b) 



a 

00 

PQ 




0.5 X 



> 



o 

oo 
O 
Q 




E - E F (eV) 



FIG. 4: (a) Minority-spin Bloch spectral functions in vac- 
uum 5.4 A above the Fe/Feo.64Nio.36 surface. The BSF is 
plotted for different fey -vectors along the path from fen = 
f = (0,0) (bottom) to fcy = 0.5 X ee (0.5,0) (top). The 
curves are shifted vertically according to their corresponding 
fc||-vectors. (b) Local density of states in vacuum 5.4 A above 
the Fe/Feo.64Nio.36 surface. 



The energy dispersion of the surface resonance at the 
iron-covered surface is illustrated in Fig. 01a). The 
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FIG. 5: Bloch spectral functions in vacuum 5.4 A above 
the surface. (a) The top surface layer composition is 
changed: i) Fe/Fe .64Ni .36, ii) Feo.75Nio.25/Feo.64Nio.36, iii) 
Feo.50Nio.50/Feo.64Nio. 36, and iv) Ni/Feo.64Nio.36. 
(b) The subsurface layer composition is changed: i) 
Fe/Fe/Fe .64Ni .36, ii) Fe/Fe .5oNio.5o/Fe .64Nio.36, and iii) 
Fe/Ni/Feo.64Nio.36- The curves are shifted vertically in or- 
der to be better distinguishable. 



surface resonance, well resolved at fcy = Y, exhibits 
strong dispersion along the T — X direction in the two- 
dimensional Brillouin zone. The LDOS is simply the 
corresponding ku -integrated BSF and it is shown in 
Fig. 0|b). The LDOS reproduces again a pronounced 
peak for a low binding energy but its position is slightly 
shifted upwards. This result follows from the fact that 
contributions from the whole surface Brillouin zone are 
included in the LDOS on equal footing. The sharp peak 
observed in the STS experiment indicates that the tip 
apex acts in this case as a sharp fc||=0-filter. The con- 
tributions from the central part of the surface Brillouin 
zone dominate. The broadening effects caused both by 
the energy dispersion of the surface resonance and the 
tip - sample interaction (not included in our analysis) are 
more important for small tip - sample distances (below 
5 A)2. 

The concentration dependence of the surface resonance 
is shown in Fig. [5ja). The strength of the surface res- 
onance is reduced with increasing content of Ni-atoms 
in the topmost surface layer: it almost disappears for 
less than 75 % Fe in the surface layer. Another surface 
resonance is formed on the Ni-rich surfaces. This reso- 
nance is shifted to lower energies (about 1 eV below the 
Fermi level) and also corresponds to d 2 2-like minority- 
spin states. We have also explored how the surface res- 
onance on the iron-rich surfaces depends on subsurface 
composition of the sample [see Fig. [3Jb)] . It turns out 
that the energy position of the surface resonance depends 
only weakly on the composition of the subsurface layer. 
We can conclude that the position of the surface reso- 
nance in the spectra is primarily determined by the com- 
position of the topmost layer and, thus, a good indicator 
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for the local concentration in this layer. This confirms 
the experimental observations. 

IV. CONCLUSION 

The chemical contrast in atomically resolved STM im- 
ages of fee Fe-Ni alloy surfaces can be explained by ab 
initio calculations which show that the (001) surface of 
Fe64Ni36 is buckled: Fc atoms are shifted outwards; Ni 
atoms inwards. We find that although the surface Fe 
concentration is w 59%, the surface is mostly c(2 x 2) 
ordered, which accommodates only 50% Fe. The excess 
Fe atoms are inserted in frequent very short Fe-rich anti- 
phase domain boundaries. 

In contrast to surface alloys based on a homogeneous 
bulk, however, the particular distribution of subsurface 
Fe and Ni atoms additionally modulates the apparent 
height of the surface atoms in the STM images. Com- 
pared to other alloys, this effect is rather strong on the 
Fee4Ni36(001) surface, suggesting significant deviations 
from a random subsurface distribution, e.g., a decompo- 
sition into regions with 50% Fc and 100% Fe similar to 
the structure of the surface monolayer. 

We have also observed a distinct peak in differential 
tunneling conductance spectra at a negative sample bias 
of roughly 0.35 eV. However, this peak is strong only in 
rather small regions on the surface. On the basis of our 
ab initio electronic structure calculations, we attribute 
this peak to a d z 2 -derived surface resonance band that is 
formed at the T-point of the surface Brillouin zone and 
in its close vicinity with its lower edge 0.3 eV below the 



Fermi level. The resonance occurs in model surfaces with 
high Fc content (75 % and more in the topmost layer) and 
therefore indicates regions on the surface which are Fe 
enriched, e.g., the mentioned Fe-rich domain boundaries 
of the c(2 x 2) ordered surface. 

These results show that with certain limitations both 
surface and subsurface chemical information can be ex- 
tracted from STM and STS images of the Fe 64 Ni 36 (001) 
surfaces: (1) Constant current images with atomic res- 
olution at very low tunneling voltages provide chemical 
contrast in an atom-by-atom fashion, which, however, is 
mixed with a height modulation corresponding to the 
subsurface composition. (2) Constant current images at 
high, in particular positive, tunneling voltages, which 
are less sensitive to the surface resonance intensity, 
map an apparent height, which is a combined effect 
of the composition of the topmost two or even more 
monolayers. (3) Current images at voltages just below 
the surface resonance peak, which are very sensitive 
to the surface resonance intensity, show the surface 
composition with a lateral resolution on the nanometer 
scale without significant interference by the subsurface 
composition. 
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